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Abstract 

We propose a three-pulse coherent ultrafast optical technique that is particularly sensitive to 
two-exciton correlations. Two Liouville-space pathways for the density matrix contribute to this 
signal which reveals double quantum coherences when displayed as a two-dimensional correlation 
plot. Two-exciton couplings spread the cross peaks along both axes, creating a characteristic highly 
resolved pattern. This level of detail is not available from conventional one-dimensional four-wave 
mixing or other two-dimensional correlation spectroscopy signals such as the photo echo, in which 
two-exciton couplings show up along a single axis and are highly congested. 

PACS numbers: 78.47.+p,78.67.De, 42.50.Md,71.35.Cc 
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Investigating the correlations of multiple excitons in semiconductors is a challenging 
many-body problem that had drawn considerable theoretical [1, 2, 3, 4, 5] and experimental 
[6, 9] attention. Correlations of two-excitons beyond Hartree-Fock (HF) may result in either 
a red-shift (two-exciton binding energyTBE) or a blue-shift (two-exciton scattering energy, 
TSE). In GaAs semiconductor quantum wells both couplings are a few meVs and may not 
be easily resolved. Two-exciton formation has been reviewed in Refs.[6, 7, 8]. Coherent 
ultrafast four-wave mixing (FWM) [6] provides a direct probe for two-exciton correlations 
in quantum wells. The best evidence for bound two-excitons and the most accurate TBE 
in GaAs quantum wells are obtained by time integrated FWM (TIFWM), where signals are 
displayed as a function of a single (time or frequency) variable. [6, 10, 11] Quantum beats 
were observed in the negative- delay two-pulse signal along — where k^ arrives first, 
and in the positive-delay three- pulse signal along ka -|- k(, — kc (kc comes after k^) . Many 
attempts have been made to enhance the spectral resolution by displaying FWM signals 
versus two time variables. [12, 13, 14] 

Multidimensional analysis of coherent signals is commonly used in NMR to study cor- 
relations between spins. [15] These techniques were recently extended to the femtosecond 
regime[16] and applied to several chemical and biological systems. [17, 18] Three ultra- 
short laser pulses generate a signal which is heterodyne-detected by a fourth pulse in a 
chosen phase-matching direction (Fig. 1 top left). Three time delays (ti,t2,i3) can be con- 
trolled between the chronologically ordered pulses, ki, k2, ks and the heterodyne pulse k^. 
For an excitonic system the signal can be generated along the phase-matching directions, 
— ki + k2 + ks, ki — k2 -I- ks and ki -|- k2 — ka. We shall denote these three techniques 
as 5*1, S\\ and S\\\ [19, 20] respectively. The time-domain nonlinear response is given by 
combinations of multi-time correlation functions which depend on the time delays ti, and 
^3. Displaying the signal as a function of two time delays (or their conjugate frequencies) 
while holding the third fixed gives the two-dimensional correlation spectroscopy (2DCS) sig- 
nals. Note that ti, ti and tz are positive. This is different from conventional FWM where 
there is no fixed time ordering among the k^, k^, k^ pulses and their delays can be either 
positive or negative. We shall show that controlling the time ordering in 2DCS is critical for 
the unambiguous identification of many-body correlations. 5*1 (photon echo) signals were 
recently reported in GaAs heterostructures[21, 22, 23, 24]. In this work, we demonstrate 
that the technique, when displayed in a proper projection plane, can access electronic 
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correlations beyond HF such as bound and unbound two-excitons with very high resolution, 
not available in ID FWM and any other 2DCS techniques. 

We employ a multi-band ID tight-binding Hamiltonian[25, 26] for calculating the 
signal. It includes a simplified Coulomb interactions and accounts for the heavy-hole (HH), 
light-hole (LH) excitons and their continuum states in a tractable way. This model has been 
successfully used to reproduce many ID FWM[26] and the recent Si 2DCS[23] experiments. 
For clarity we focus on the HH spectral region. With near-resonant excitations, only 
two terms contribute to the third-order optical response in the phase-matching direction 
ki -I- k2 — ks. These are represented by the double-sided Feynman diagrams shown in the 
bottom of Fig. l(lcft), where e and / represent respectively HH exciton and two-exciton 
manifolds and g represents ground state (top right of Fig. 1) The rules of these diagrams 
are given in Ref . [19] . 

During ti and t2, both pathways are respectively in a single (peg) and double (pfg) quan- 
tum coherences. They differ only in the third interval ts {peg for (i) and p/e for (ii)). 
Signatures of two-exciton couplings thus show up along ^3 and t2 but not ti. To enhance 
the resolution, we therefore look at (t3,t2) correlation plots. We shall display the signal 
by a double Fourier transform with respect to ^3 and t2 with the corresponding conjugate 
frequencies fl^ and fl2, holding ti fixed. The expected 2DCS, Sm {fl3,fl2,ti), is schemat- 
ically sketched at the bottom of Fig. 1. The individual contributions of pathways (i) 
and (ii) are denoted by the symbols given at the bottom of each Feynman diagram. / 
can be either red-shifted (bound) / = 2e — A^, un-shifted f — 2e or blue-shifted (scat- 
tering) f = 2e + relative to twice the HH exciton energy, 2e. Red(blue) energy shift 
corresponds to two correlated excitons with opposite (same) spins. Open circles represent 
un-shifted two-excitons and solid (open) symbols represent red(blue) shifted two-excitons. 
Energy shifts from dressed states are negligible in the present weak-field x^^"* limit and thus 
the HF (mean-field) contribution gives energy-unshifted two-excitons. 

We now examine the five peaks in Fig. 1. We note that the peaks resulting from (i), 
are only spread along ^2 and form a vertical pattern since it has eg resonances along ^3, 
while for (ii) the peaks are spread along both axes (diagonal pattern). Four properties 
make this signal particularly useful for resolving various two-excitons and even the detailed 
structure of the two-exciton continuum. First, along t2 there are only fg resonances and 
thus the 0,2 axis provides a clean projection for two-excitons, without interference from 
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eg resonances. Second, the two-exciton correlation energy can be obtained from a single 
peak. is extracted from the solid square at (e, 2e — A J and A^ from the open square 
at (e,2e + A J. In TIFWM, in contrast, A^^ is determined by subtracting the energies of 
two peaks. However, the HF peak (the circle) is not always clearly visible for certain pulse 
polarizations and laser detunings. Third, the squares from (i) for different / are generally 
much stronger than triangles from (ii). Thus it is sufficient to get the two-exciton energies 
from the well-resolved and strong (i) peaks (squares) along 0,2 even when the triangles are 
not resolved. Fourth, both axes ft2 and O^ involve two-exciton resonances. Along Q2, 
there are only fg resonances, while along Q3 there is both eg (i) and fe (ii). Spreading 
two-exciton resonances along both axes is critical for the high resolution of two-excitons. 
Along the two-exciton contribution (/e in (ii) ) is weaker than the single-exciton one {eg 
in (i)) and thus may not be easily resolved. However, we shall show later that even this 
weak signature along Q3 is crucial for achieving high resolution when combined with the 
information from the 0,2 projection. 

Using the tight-binding Hamiltonian[25, 26] with 10 sites, we have implemented the x*^^^ 
formalism for closing the infinite hierarchy of equations that has been widely applied to 
Frenkel[27, 28] and Wannier[2, 29] excitons. In Fig. 2, we present the simulated 2D spectra 
obtained by solving Eqs. 20 and 21 of Ref.[22]) using periodic boundary conditions. Using 
the same parameters of Ref. [30], we set the dephasing times for excitons Tex = 2 ps and 
two-excitons T2ex = 1 ps in panels A-E. Panel E is obtained by cross-circular Gaussian 
pulses, where is hnearly polarized and ks, k2 and ki are respectively right, left and right 
circularly polarized. All other panels are calculated with co-hnear Gaussion pulses. Pulse 
parameters are chosen such that red shifted two-excitons are selectively excited. The optical 
power spectra width, Aopt, and carrier frequency, u;c(relative to e) are given in the caption. 

In all panels, the origin is e for ^3 and 2e for O2. All peaks or shoulders are assigned 
using the symbols given in Fig. 1. Circle and solid square are strong and well resolved along 
O2, while sohd square and triangle are not resolved along O^. The O2 value of solid square 
gives the TBE, A^. Panels A and B are calculated respectively with and without correlated 
two-excitons. As expected, features related to red shifted two-excitons such as the solid 
square and triangle, disappear. Panels C and D repeat the calculations of A and B for 
pulses tuned to the blue-shifted two-excitons. One can obtain blue-shifted TSE, A^ , from 
the open triangle (a shoulder but not a resolved peak) in panel C, even though we can not 
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resolve open square from two-exciton continuum. In panel E (cross-circular excitation) , the 
HF contribution (the circle) at the origin (e, 2e)=(0,0) is missing. However, one can still 
obtain the TBE using the sohd square. We also note that TBE may be extracted despite 
the large line broadening. Panel F shows the 2D spectrum calculated with faster dcphasing 
times Tex = 1-0 ps and T2ex = 0.5 ps. Even though we can no longer resolve circle and solid 
square along the clean axis, it is still possible to obtain their energies from the 'ridge' of 
the contour lines (two white lines). The TBE is given by the fl2 value of the bottom hue. 

The TBE can be obtained from the intense and well- resolved solid square in panel A. We 
thus circumvent the difficulty of probing TBE from the splitting between solid square and 
the unresolved solid triangle, which overlap with eg and fe resonances along ^3. Additional 
calculations also show that the TBE related solid square always has a solid triangle to the 
red, and the TSE related open square always has a open triangle to the blue, as predicted in 
Fig. 1. Moreover, the solid square and triangle (open square and triangle) always appear 
in pairs and thus we can easily obtain accurate TBE(TSE) by identifying both features, 
even when line broadening is large (e.g. panel F). The much weaker solid triangle from 
pathway (ii) plays a crucial role in identifying the solid square in panel F. Similarly, in panel 
C we can not identify the open square without the help of the open triangle. Thus it is 
impossible to resolve the signature of blue-shifted two-excitons (unresolved open square and 
triangle) in panel C with conventional ID FWM where the signature will either be covered 
by the two-exciton continuum along Q2 or along Q3 by the broadening of the much stronger 
single-exciton peak, the circle. This signature can not be obtained from any other 2DCS 
techniques where two-excitons show up along a single axis. However, with the panoramic 
2D view offered by Sm, we can easily obtain very accurate TSE, even though the open 
square and triangle are weak may not be resolved along any single axis. 

In Si, we had demonstrated the first partial separation of two-excitons[22] along O3 due 
to the overlapping eg and fe resonances along ^3, as shown by the elongation of peaks in both 
experiments and simulations. [22, 23] Sm provides an additional separation of two-excitons 
by spreading them along ^2- Without this separation, one can not resolve any two-exciton 
features in panel F, let alone the accurate TBE. It is the combination of the two dimensions 
fl^ and fl2 that makes it possible to go from two ambiguous, unresolved open square and 
triangle in panel C, to retrieve the unique signature of blue-shifted two-excitons and obtain 
the TSE. 
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The Siii technique provides a new perspective into the capacity of TIFWM experiments 
to provide TBE.[10, 11] Let us recast the existing TIFWM signal using 2DCS terminology. 
The two-pulse signals[10] along 2kb — is given by Wa (h) — \ Si (^3,^2 = 0,ti)f dt^ 
(positive delay) and Wb (^2) = \Sui (isjhjti = 0)f dt^ (negative delay). Three-pulse 
TIFWM signals[ll] along k„ + k;, - k^ correspond to Wc (^2) = Jl^ \Sui (^3, ^2, = t?) |^ dts 
(positive delay) and Wd (^2) — \Su (^3, h-, h — i?)|^ dts (negative delay). Strong quan- 
tum beats only show up in Wb and Wc, both related to and depend on t2. This is 
clear from our pathway analysis. First, along t2 there are only two-exciton resonances (fg) 
as shown in diagrams (i) and (ii) of Fig. 1, and thus we have well-defined quantum beats 
along ^2 and well-resolved peaks along Second, the squares from (i) show up mainly as 
eg resonances along ^3 and can be much stronger than the triangles, although they are all 
induced by correlated two-excitons. Without such correlations, there will only be a circle 
at the origin. It is the stronger solid square and circle that give quantum beats. However, 
there are no appreciable beats for Wa and Wd- Fig. 3 shows the pathways for Sj and the 
schematic 2D spectra projected respectively onto (Q3, —Qi) and 0.2)- Obviously, there 
are no two-exciton splitting along either fli or ^2 and thus Si does not show quantum beats 
with respect to ti or t2, as in Wa- A similar conclusion applies to 5*11 , which corresponds 
to Wd- The only axis along which quantum beats occur for and is ^3. However, it 
is hard to observe such beats with any technique (including Sm), because there are always 
dominant eg resonances along ^3 for all three 2DCS techniques and thus the two-exciton 
resonances are not well resolved. 

Siii also demonstrates an important limitation of TIFWM in probing TBE. Our calcu- 
lations show that TIFWM only yields the correct TBE in ideal cases such as Fig. 2A where 
the ^^2 coordinate of the circle coincides with 2e and its intensity is comparable to solid 
square. However, the ^2 value and intensity of the circle are very sensitive to the pulse 
width, detunings, pulse polarizations, and are further affected by the cancellation of two- 
exciton continuum to the HF portion. Calculations with long dephasing times (not shown) 
also suggest that the two-exciton continuum is far from uniform. Although the solid square 
has a fixed Q2, its intensity is also very sensitive to the pulse properties. Therefore the circle 
and solid square do not always dominate the signal. Other components from two-exciton 
continuum or even open square may also contribute or even dominate the quantum beats. 
This makes the beating experiment less accurate for probing TBE. Furthermore, in the 
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most favorable case for probing TBE where only bound two-excitons are generated (e.g. the 
cross-circular excitations in Fig. 2E), the HF peak (circle) disappears completely. In this 
case, any beating frequency obtained is not connected to TBE. Therefore, one has to use 
other pulse polarizations to obtain quantum beats where the position and intensity of the 
circle is affected by other types of two-excitons. 

The advantage of Sm is even more pronounced when HH and LH excitons both contribute 
and TIFWM signal is congested with many types of beating components along t2 arising 
from pure HH (LH) and mixed two-excitons. By dispersing all pathways into an extra 
dimension and spreading them along Q2 over a spectral range twice that of Si and Sn 
techniques, Sm achieves a complete separation of pure HH (LH) two-excitons from the mixed 
two-excitons. Furthermore, by choosing specific pulse polarizations, the contribution from 
correlated pure or mixed two-excitons can be completely isolated from their corresponding 
HF contributions. 

In summary, we propose a 2DCS technique, Sm, that is intrinsically sensitive to two- 
exciton correlations and can distinguish between different species of excitons and two- 
excitons in photo-excited semiconductors by spreading them along two axes. Combining the 
information from different Liouville space pathways, [19] Sm can accurately retrieve TBE 
even when it is smaller than exciton line broadening. We further report an unambiguous 
signature of correlated unbound two-excitons lying underneath the two-exciton continuum. 
Conventional ID or 2D techniques can not achieve such resolution for the lack of time order- 
ing among pulses, which is essential for the identification of specific pathways. For example, 
either ID or 2D pump- probe signals generally mix the 5*1 and Su techniques and thus one 
can not employ them to identify different pathways, let alone combine them to give higher 
resolution. Sm also provides new insights into the quantum beats in the most accurate 
TIFWM for obtaining TBE[10, 11] of quantum wells, which are dominated by one of the 
two pathways and generally do not provide accurate TBE. 

Sm could be applied to other strongly correlated systems such as Bose- Einstein conden- 
sate and superconductors. For instance, it could be helpful in the considerable current 
effort aimed at the coherent manipulation of excitons and exciton Bose-Einstein condensa- 
tion. Two-exciton effects arc critical for achieving these goals. 

The support of the Chemical Sciences, Geosciences, and Biosciences Division, Office of 
Basic Energy Sciences, U. S. Department of Energy is gratefully acknowledged. 
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Figure captions 

Fig. 1 Top: Schematic experimental setup (left) and exciton level scheme (right). Bot- 
tom: Feynman diagrams of Sui technique (left) and schematic 2D spectrum (right). 

Fig. 2 (Color online) Calculated 8111(^23,^22,^1) for (A) and (B), Aopt = 3.9 meV, 
Uc = -2.85 meV; (C) and (D), Aopt = 1-8 meV, Uc = 2.15 meV; (E), Aopt = 3.9 meV, 

= -3.85 meV; and (F) Aopt = 3.9 meV, = -1.85 meV. 

Fig. 3 (Color onhne) Top: Feynman diagrams of 5*1 technique. Bottom: Schematic 2D 

spectra, Si{Q:s,t2, -^1) (left), Si{^l3,Q.2,ti) (right). 
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